function y = CL_deltaee(x)

  if x >= 0
      c = [0 0 5.0064 -6.1541 1.5511 0.4677 0.0002];
      p = [x.^ 6, x.^ 5, x.^ 4, x.^ 3, x.^ 2, x.^ 1, x.^ 0]';
  else
      c = [0 0 -5.064 -6.1541 -1.5511 0.4677 -0.0002];
      p = [x.^ 6, x.^ 5, x.^ 4, x.^ 3, x.^ 2, x.^ 1, x.^ 0]';
  end
  
  y = c * p;
  
end
